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An  algorithm  for  solving  the  problem  of  finding  a 
maximal  dynamic  flow  through  a  network  Is  described.  No 
proofs  are  given. 
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CONSTRUCTION  OF  MAXIMAL  DYNAMIC  FLOWS  IN  NETWORKS 
L.  R.  Fordf  Jr.«  and  D.  R.  Fulkerson 

1.  Introduction.  In  this  note  we  describe  briefly  an 
easy  and  efficient  algorithm  for  solving  the  following 
problem.  Suppose  given  a  network  (linear  graph)  In  which 
each  link  has  associated  with  It  two  positive  Integers,  one 
a  coQBiodlty  flow  capacity,  the  other  a  traversal  time. 

Assuming  that  some  node  of  the  network  Is  a  source  for  the 
comaodlty,  another  a  sink,  and  the  remaining  may  either 
transship  the  coniiiodlty  Immediately  on  receipt  or  hold  for 
later  shipment,  what  Is  the  maximal  amount  that  can  be  shipped 
from  source  to  sink  In  any  given  number  of  time  periods? 

For  example.  In  the  toy  network  of  Figure  A,  P^  Is  the 
source,  P^  the  sink,  and  the  capacities  of  the  links  In  terms 
of  flow  per  unit  time  are  the  first  numbers  of  the  pairs,  the 
traversal  tlises  the  second  numbers.  How  many  units  of 
comodlty  flow  can  reach  P^  from  Pq  In  five  time  periods, 
say,  and  what  Is  a  flow  pattern  which  achieves  this? 
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A  simpler  problem  related  to  the  dynamic  one  just 
described  is  the  maximal  static  flow  problem,  i.e.,  assuming 
a  steady  state  condition,  find  a  maximal  flow  from  source  to 
sink  in  a  capacitated  network.  The  labeling  process  [7] 
provides  a  simple  method  of  solution  for  this  latter  problem; 
a  variation  of  this  process  is  also  used  as  a  subroutine  in 
the  algorithm  for  the  dynamic  problem. 

The  construction  of  an  optimal  dynamic  flow  for  a  given 
number  T  of  time  periods  first  produces  a  static  flow  in  the 
network  which  has  certain  prescribed  properties.  The  dynamic 
solution  is  then  obtained  from  the  static  solution  by 
decomposing  the  latter  into  "chain  flows,"  starting  each 
chain  flow  at  time  zero,  and  continuing  each  so  long  as  there 
is  enough  time  left  in  the  T  periods  for  the  flow  along  the 
chain  to  arrive  at  the  sink.  Thus,  for  the  network  of 
Figure  A,  the  static  flow  produced  by  the  algorithm  for  T  -  5 
is  shown  in  Figure  B  below.  This  flow  decoBq;>08es  into  a  flow 
of  two  units  along  each  of  the  chains  ^0^2^3‘ 

each  of  these  chains  has  a  total  traversal  time  of  four,  each 
can  be  used  twice  in  the  five  periods,  hence  eight  units  of 
the  commodity  reach  from  Pq  in  the  total  time  interval  in 
question.  The  fact  that  a  maximal  dynamic  flow  of  this 
particularly  simple  kind  exists  for  all  T  seems  to  us  rather 


remarkable. 
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Plgure  B 


No  proofs  will  be  given  in  this  paper.  These  may  be 
found  in  [9] . 

2.  Definitions  and  Heuristic  Discussion.  Let 

Tq,  T^0  the  nodes  of  a  network  N,  and  denote  by 

PiPj  the  link  from  to  Fj  (in  that  order).  Associated  with 

each  link  P^Fj  present  in  N  are  two  positive  integers, 

and  t^j,  its  capacity  and  traversal  time,  respectively.  (In 

the  network  of  Figure  A  the  capacities  and  times  are  assumed 

aynmetric,  e^j  -  c^^  and  t^^  -  t^^,  but  this  is  not  necessary 

in  the  algorithm.)  The  notation  is  taken  so  that  Fq  is  the 

source,  F^  the  sink, 
n 

A  static  flow  in  N  is  a  collection  of  numbers  (x^j)« 
one  for  each  link  F^Pj  of  N,  such  that 

(1)  2  (x..  -  X..)  -  0  (i  =  1#  n-1) 

J-0 


°  ^  *ij  ^  'ij  • 
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(It  is  assumed  in  (1)#  and  elsewhere^  that  a  variable 
appears  only  if  ^  link  of  N.)  The  value  of  the 

static  flow  (x^j)  is  the  total  amount  leaving  Pq# 

Z  (Xqj  •>  Xjq)  .  Thus  the  static  maximal  flow  problem  may  be 


J 


expressed  as 


(2) 


n 


maximize 


subject  to  the  equations  and  inequalities  (1).  In  this  form# 
I  the  static  problem  can  be  solved  by  the  algorithm  given  in 
[7]. 


An  alternative  formulation  of  the  static  problem  appears 
if  one  focuses  attention  on  chains  in  N  which  lead  from  Pq 
to  P^  [6].  Thus,  let  C^,  be  a  listing  of  chains  from 

r©  to  p„.  ....  Lp  a  listing  of  links  of  N,  and  suppose 

^1'  *'*'  \  represent  the  amounts  of  flow  along  the  chains. 

If  (Spg)  is  the  p  by  m  incidence  matrix  of  links  vs.  chains, 
then  the  problem  is  to 

m 


(3)  maximize  Z  x. 

k-1  ^ 


subject  to  the  restrictions 


where  c^  is  the  capacity  of  the  r-th  link. 


The  chain  foxmilation  of  the  static  problem  does  not 
seem  to  be  fruitf\)l  in  a  computational  sense  directly. 

However#  looking  at  the  problem  in  this  form  provides  sosie 
insight  into  the  dynamic  problem.  Suppose  one  knew  a  priori 
that  there  were  a  maximal  dynamic  flow^  of  the  kind  described 
in  the  introduction.  Let  be  the  traversal  time  of  the 
cltain  C^.  Then  the  problem  would  be  simply: 

m 

(5)  maximize  2  x.  (T  +  1  -  t^) 

k-1  ^  ^ 

subject  to  the  constraints  (4)«  since  for  ^  T  +  1#  the 
coefficient  T  1  -  t^  is  the  number  of  times  can  be  used 

in  the  T  periods#  and  for  tj^  >  T  +  1,  a  solution  which 
maximizes  assigns  x^  0. 

How  problem  (4)  —  (5)  can  be  put  back  in  a  form  similar 
to  that  of  (l)-(2)#  i.e.# 

(6)  maximize  (T+l)  (*oj  "*J0^  ”  ^  ^Ij^iJ 

J  1#  J 

subject  to  (l).  It  is  not  difficult  to  see  that  this  may  be 
restated  as  a  capacitated  Hitchcock  transportation  problem# 
for  which  several  single  computational  methods  have  been 

proposed  [l#  2#  8J . 

To  recapitulate#  if  one  could  establish  the  existence  of 
a  maximal  dynamic  flow  of  the  simple  kind  previously  described# 
the  dynamic  problem  would  not  be  difficult  computationally. 

^e  shall  give  no  formal  definition  of  a  dynamic  flow  in 
N#  since  the  algorithm  of  the  paper  deals  essentially  only 
with  static  flows. 
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Encouraged  by  the  fact  that  one  of  the  authors  (L.  R.  Ford*  Jr.) 
had  constructed  the  existence  proof  needed  on  the  assun^tlon 
that  T  was  sufficiently  large,  we  were  led  to  an  examination 
of  the  primal-dual  method  [3,  6]  for  solving  the  problem  of 
finding  a  maximal  static  flow  which  minimizes  total  flow  time 

Z  t_x. .  over  the  class  of  all  maximal  static  flows,  i.e^, 

1,J 

of  solving  (1)  and  (6)  for  large  T.  Surprisingly,  it  turned 
out  that  the  algorithm  was  grinding  out,  sequentially,  maximal 
dynamic  flows  for  T  •  0,  1,  2,  ...,  and  that  the  proof  of 
this  fact  could  be  given  directly  from  the  algorithm  itself. 


3.  Dynamic  Algorithm.  Routine  I  of  the  algorithm  is 
an  iterative  process  which  constructs  an  integral  static 
flow  together  with  nonnegative  integers  one  for 

each  node  having  the  properties 


(7)  a.  Tq  -  0,  -  T  +  1 

b.  -> 

’"l  ''  ’■j  *1J 


0 

""ij 


d. 


To  state  the  routine,  we  suppose  we  have  an  integral 
flow  (x^j)  and  node  Integers  satisfying  (7)  with  “  t, 
and  wish  to  construct  (x^j)  and  satisfying  (7)  with 

I 

»  t  +  1.  To  start  out,  one  may  take  all  x^j  »  0,  all 
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Llnks  for  which  -  Tj  are  called  admlaalble 

below.  Notice  that  at  most  one  of  P^Pji  PjPj^  will  be 
adinl88lbl6«  and  that  Initially  no  llnka  are  adnlaalble. 


h 


Routine  I. 

(a)  Label  Pq  with  (P^«  oo);  conalder  Pq  aa  unscanned. 

(b)  Take  any  labeled,  unscanned  node  P^  (initially  Pq  will 
be  the  only  such  node);  suppose  labeled  h).  To 


all  nodes  P^  which  are  unlabeled  and  such  that  P^P^ 

Is  admissible  and  assign  the  label 

(P^,  min  (h,  c^j  -  To  all  nodes  Pj  which  are 

now  unlabeled  and  such  that  P^P^  Is  admissible  and 
Xji  >  0,  assign  the  label  (P^,  min  (h,  Consider 

P^  as  scanned  and  newly  labeled  P^  (If  any)  as 
unscanned.  Repeat  until  P^  Is  labeled  or  until  no 


new  labels  are  possible  and  P^  Is  unlabeled.  In  the 
former  case  go  to  (c)  below;  In  the  latter  case,  let 
(x^j)  denote  the  flow,  and  proceed  to  (d). 

(c)  If  Pj^  Is  labeled  (Pj^,  h),  replace  x^^  by  +  h;  If 
P^  is  labeled  (P",  h) ,  replace  x^  by  -  h.  In 


either  case,  next  turn  attention  to  Pj^.  In  general. 
If  Pj^  Is  labeled  (Pj,  in),  replace  x^^  by  Xjj^  +  h,  and 


If  labeled  (Pj,  ra),  replace  x^^j  by  x^^j  -  h.  In  either 


^There  are  more  efficient  ways  to  start  the  routine,  e.g., 
the  procedure  of  [5]  for  finding  a  "shortest  chain  from  Pq  to 
may  be  used. 


n 
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case  tuming  attention  then  to  Pj,  Stop  the  replace¬ 
ment  when  Pq  is  reached.^  Starting  with  the  new 
integral  flow  thus  generated »  discard  the  old  labels 
and  repeat  (a)  and  (b)  until  the  latter  case  of  (b) 
obtains . 

I 

Define  by 


if  is  labeled* 

+  1  if  Pj^  is  unlabeled. 


Repeat  routine  I  starting  with  and  (giving  as 

new  admissible  links  those  P^Pj  for  which  +  t^j  *  Xj)  until 

the  value  of  has  been  boosted  to  T  -i-  1.  At  this  point  a 
solution  of  the  dynamic  problem  for  T  periods  is  at  hand 
merely  by  decomposing  the  final  flow  obtained  into  chain 
flows.  Routine  II  below  uses  a  labeling  procedure  to  effect 
this  decomposition. 


Routine  II. 

(a)  Label  Pq  with  (P^,  oo ) ;  consider  Pq  as  unscanned. 

(b)  Take  any  labeled,  unscanned  node  P^^;  tuppose  labeled 
(P,.,  h).  To  all  nodes  P.  which  are  unlabeled  and  such 

iC  j 

that  >  0,  assign  the  label  (P^,  min  (h*  x^^)). 
Consider  as  scanned  and  newly  labeled  Pj  (if  any) 
as  unscanned.  Repeat  until  either  P^  is  labeled  or 

^The  static  flow  value  has  been  increased  by  h  at  this 
point. 
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new  labels  are  iapoasible  and  Is  unlabeled.  In  the 

former  case,  proceed  to  (c);  in  the  latter  case,  stop. 

(c)  If  Pjj  Is  labeled  (Pj^,  h),  replace  by  -  h,  next 

turning  attention  to  Pj^.  In  general.  If  Pj^  Is  labeled 

(Pj#  n),  replace  Xjj^  by  -  h,  and  proceed  to  fy 

4 

Stop  the  replacement  when  Pq  is  reached. 

Repeat  routine  II  unti3  the  latter  case  of  (b)  obtains; 
the  flow  (x^j)  has  then  been  decomposed  into  chain  flows. ^ 

The  decomposition  obtained  in  this  way  is  of  course  not  unique 
in  general.  However,  any  list  of  chain  flows  obtained  from 
the  final  flow  of  routine  I  represents  a  solution  to  the 
dynamic  problem  for  T  periods.  If  the  k~th  chain  of  the  list 
has  a  flow  of  x^^  and  a  traversal  time  of  t^^,  then  the  maximal 
dynamic  flow  value  is  given  by 

Z  Xj^(T  +  1  -  tj^)  . 

Ic 

If  one  is  interested  only  in  this  value  and  not  the  flow 
pattern,  routine  II  of  the  algorithm  may  be  ignored.  Instead, 
let  be  the  final  node  numbers  produced  (l.e.,  -T  1), 
and  define 

(8)  -  max  (0,  Tj  -  -  t^j)  . 

4 

One  sees  that  at  the  end  of  process  (c)  a  chain  carrying 
h  units  of  flow  has  been  traced  out  (in  reverse). 

^or  an  arbitrary  flow,  circulations  and  chain  flows  from 
Pjj  to  pQ  nay  remain.  'Riis  cannot  happen,  however,  for  a  flow 

produced  by  routine  I. 
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Then  It  can  be  shown  that  the  maximal  dynamic  flow  value  is 
given  by 


(9) 


The  node  numbers  ir^  and  link  numbers  provide  other 
significant  information  also.  For  example,  links  with 
positive  are  "bottlenecks"  for  the  dynamic  flow  —  more 
specifically,  a  link  P^Pj  for  which  Yj^j  ^  0  bottlenecks  the 
flow  in  time  periods  1,  2,  . . . ,  ir^  +  ^ij’^  Thus  the 

flow  value  for  T  periods  can  not  be  increased  unless  some  one 
of  these  links  is  given  a  larger  capacity. 

Notice  that  the  dynamic  algorithm  has  the  feature  alluded 
to  in  the  introduction,  namely,  that  in  solving  a  problem  for 
a  given  number  of  time  periods  Tq,  optimal  solutions  for  all 
fewer  numbers  of  time  periods  T  <  Tq  can  be  obtained  as  by¬ 
products.  In  addition,  it  can  be  shown  that  at  some  stage  T 
in  the  computation,  a  general  solution  of  the  problem  for  all 
numbers  of  time  periods  T  ^  T  is  obtained.^  (The  value  T  is 
never  greater  than  the  maximal  traversal  time  from  Pq  to  P^^, 
and  is  usually  much  less.)  At  the  stage  T,  the  static  flow 
obtained  is  a  maximal  flow  which  minimizes  total  flow  time 

Z  t.  .X. .  over  all  maximal  static  flows,  and  the  links  which 
i,J 

lead  from  the  final  labeled  set  of  nodes  to  its  complement, 

^See  the  definition  of  the  minimal  dynamic  cut  in  [9]  . 

7 

'Thus  the  solution  of  the  problem  for  all  T  is  a  finite 
process. 
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the  unlabeled  8et«  are  a  miniraal  (atatic)  cut  [4,  6,  7]  in 
the  network. 
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